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Abstract. We consider the relationship between Hele-Shaw evolution with drift, the porous medium 
equation with drift, and a congested crowd motion model originally proposed by [MRS]- [MRS V] . We 
first use viscosity solutions to show that the porous medium equation solutions converge to the Hele- 
Shaw solution as m — > oo provided the drift potential is strictly subharmonic. Next, using of the 
gradient flow structure of both the porous medium equation and the crowd motion model, we prove 
that the porous medium equation solutions also converge to the congested crowd motion as m — > oo. 
Combining these results lets us deduce that in the case where the initial data to the crowd motion model 
is given by a patch, or characteristic function, the solution evolves as a patch that is the unique solution 
to the Hele-Shaw problem. While proving our main results we also obtain a comparison principle for 
solutions to the minimizing movement scheme based on the Wasserstein metric, of independent interest. 
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Oh 1. Introduction 

< 

Let f^o be a compact set in M. d with locally Lipschitz boundary, and let $(a;) : M. d — > E be a C 2 
"ti function which satisfies 

£h (Al) A$ > iii R d . 

For J7o and $ as given above, we consider a function u : M. d — > R, u(x, t) > solving the following free 
t-H boundary problem: 

> ( 

(N -A«(.,t) = A* in {u>0}; 

Q> I V = -d v u - d v <& on d{u > 0}. 

Here v Xtt is the outward normal vector of the set VL t {u) :— {x : u(x,t) > 0} at x € T t (u) := dflt(u), and 
V denotes the outward normal velocity of Qt{u) at x € T t (u). 

In terms of u, v = — Vu/|Vu| and thus one can write down the second condition of (P) as 
L* ut= |Vu| 2 + Vu ■ V$ on d{u > 0}, 



X 



given that |Vu| ^ at the boundary point. Note that the free boundary velocity V may be positive or 
negative depending on the behavior of $ on r(t). Consequently Qt(u) n^ay expand or shrink over time 
(see Figure 1). Indeed formal calculations based on (P) yield that Sl t preserves its volume over time. 
The initial data u(x, 0) = uq is the unique function satisfying 

— Auq = A$ in the interior of Qq, uq — on Qq. (1-1) 



Note that, due to (Al), uq is positive in flo and thus (1.1) is well-defined. Still, even starting from a 
smooth domain iloj the solution of (P) can develop finite-time singularities as its support goes through 
topological changes such as pinching and merging, and thus it is necessary to consider a notion of weak 
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v = -d v u - a,$ 




equilibrium profile 
s ( x ) = (C- $(z))_ 



Figure 1. Evolution of the positive phase, converging toward the equilibrium 



solutions. We will use the notion of viscosity solutions for (P), see section 2 for definitions and properties 
of u. Let us mention that the usual variational inequality formulation for weak solutions of Hele-Shaw 
flow, introduced by IE Jj , does not apply here due to the non-monotonicity of solutions in time variable. 

In the context of fluid dynamics, the problem (P) describes a flow in porous media. Indeed if we 
denote by u — u{x, t) the density of a fluid and define the velocity of the fluid as 



U 



_V$ - Vm, 



where V$ is the external velocity field given by $, then (1.2 1 and the incompressibility condition 

V-E/ = 
yields (P). 



(1.2) 



(1.3) 



When $ = and there is a fixed boundary in the positive phase through which the fluid is injected, 



(1.2) and dl .3b yield the classical one-phase Hele-Shaw problem HS . In this article however, our goal 



is to derive (P) from a model problem in crowd motion with hard congestion, as described below. 



1.1. A model in congested crowd motion. Let us recall the transport problem with density con- 
straint, introduced in |MRS| - |MRS V| . Formally the problem can be written as the following: we look 
for a solution p : M. d x [0, oo) — > [0, oo) satisfying 



p t + V • (/oV$) = if p < 1, and p < 1 for all times. 



(1.4) 



The density constraint is natural in many settings, and it describes motion of congested individuals. 
We refer to the articles MRS , M RSVflSan for applications and mathematical formulations of the problem 
(1.4). More rigorously, the problem can be written as 



Pt + V • (pu) = 0, u = P Cp V$, 



(1.5) 



where P is the projection operator and C„ is the space of inadmissible velocity fields which do not 



increase p on the saturated zone {p — 1}. We refer to MRSV for further description of C p . Due to the 



low regularity of the velocity field u and the non-continuous dependence of the operator Pc with respect 



to p, classical methods to study transport equations do not apply to (1.5). Indeed the uniqueness of 
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PME-D 



r^eote^ 



-5° *% h- 



HS-D 




equal a.e. 

Theorem 1.2 



Poo 



Figure 2. This diagram is a summary of the results of Theorems |1.1| and |1.2| Here 
poo denotes the gradient flow solution in the continuum limit, which in particular is a 



solution of (1.5) 



solutions for (1.5) is an open question, and is probably false without further conditions on the solution, 



given its hyperbolic nature. 



In MRSV| , the authors study the connection between the PDE (1.5) with poo, which is the gradient 



flow of the following functional Poo with respect to the 2-Wasserstein distance: 



P>o[p] 



I R d p{x)^{x)dx 



-00 



for 
for 



< 1 
> 1. 



(1.6) 



Further, they prove that when <!> is A-convex, the gradient flow solution poo is a weak solution for (1.5) 



However, the full characterization of the solution and further qualitative properties of the solution remain 
open due to the lack of available methods to study ( 1.5 ). The connection between p x and (P) has been 



hinted, but only formally in the context of particle velocity. 



o Our contribution: In this article, our main focus is on establishing the connection between the free 
boundary problem (P) and the gradient flow of Poo in the setting of patches, i.e. when the initial data 
is given as a characteristic function of a compact set £Iq, which we denote by xn - Note that since $ is 
assumed to have a positive Laplacian, solutions tend to aggregate and thus we expect that the gradient 
flow Poo(',t) will stay as a characteristic function at all times t > 0. 

We show that the preservation of patches is indeed the case, and moreover the gradient flow solution 
poo(',t) indeed coincides with the characteristic function of the set il t , which evolves according to our 



problem (P) with the initial support Qq (see Theorem 1.2 below). This result enables us to characterize 
the evolution of poo in a unique way and also helps to understand the geometric behavior of poo- A 
summary of our results is shown in Figure [2] 

In our analysis, the main challeng e is the low regularity of poo, since a priori we only know that it is 
in CV([0,T],p2(R d )) (see Theorem 4.1V)). Thus it is rather difficult to directly study the geometric 
property of poo • Instead of trying to directly show the link between the free boundary problem (P) with 
the gradient flow p^ , we use an approximation with degenerate diffusion. It has been formally suggested 
in MRS and San that one could consider approximating the gradient flow of Poo by the unconstrained 



gradient flow problem with the energy 



E m [p] 



in 



p$ dx. 



(1.7) 
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It is well known (for example, see [O]) that the gradient flow p m associated with E m solves the porous 
medium equation with drift 

p t -V-(V(p m ) + pV<P) = 0. (1.8) 



Let us denote p m as the viscosity solution to ( 1.8 1 with initial data xn ■ We will prove that as m — > oo, 
p m on the one hand converges to xti t locally uniformly, and on the other hand converges to poo(-,t) in 
2-Wasserstein distance. Thus it follows that xn t & n d Poo niust be equal to each other almost everywhere. 
The main ingredients of the proof consist of stability results from viscosity solution theory and optimal 
transport theory, both of which rely strongly on the convexity- type conditions on $. We also obtain 
comparison results and qualitative rates of convergences; see section 1.2 for precise statements. 



1.2. Summary of results. We are now ready to state our main results. The relevant assumptions, 
besides (Al) in the introduction, are stated in the beginning of section 4. 

Theorem 1.1. Let Hq be a compact set in M. d with locally Lipschitz boundary, and consider the initial 



data uq as given in (1.1). Then the following holds: 



(a) (Theorem 3.5) Assuming (Al), there exists a unique family of compact sets fit in M. starting 



from fio such that any viscosity solution u of (P) satisfies {u(-,t) > 0} = Q t f or oil t > 0. Fur- 
thermore, let p m denote the viscosity solution to (1.8 1 with initial data xn - Then as m — > oo, 
p m converges to p :— xn t locally uniformly in M. d — dVtt at each time t > 0. 



(b) (Theorem J^.2) Assume (A2) and (A3'), and consider po £ "p2(M. ) with ||po ||oo < 1 and 
J po<&dx < M. Let p m {x,t) denote the viscosity solution of (1.8) with initial data po. Then 
there exists p^ £ Cw([0,T]','P2(^- d )) such that for any T > 0, as m — > oo, p m (-,t) converges to 
Poo{'i~t) in 2-Wasserstein distance, uniformly in t for t £ [0, T], with the following convergence 
rate: 

C{M,T, || ASH*,) 



SUp W 2 {p m {t),Poo{t)) < 

te[o,T] 



,1/24 



Combining Theorem |1.1| (a) and (b), we immediately draw the following conclusion for the identifica- 
tion of Poo . 

Theorem 1.2. [Characterization of poo] Let 0,$, p^ and p as given in Theorem \l.l\ If (Al), (A2) 
and (A3') hold and if po = Xn o; then p^ — p a.e. 



As a by-product of our analysis, we also show that a version of comparison principle holds between 
solutions to the discrete Jordan-Kinderlehrer-Otto (JKO) steepest descent scheme: 

Theorem 1.3 (Comparison principle, see Theorem |5.1[ ). Let $ satisfy (A3). For 2 < m < oo, consider 
the two densities poi S 'P2,M 1 (^- d ), P02 S T J 2,M 2 (^' d ) (Pi,Mi * s a $ defined in section 5.1) with the property 



Mi < M2 and poi < P02 o,.e. (In the case m = 00, we require in addition that ||/?o* || 00 < 1 f or i — 1- 2J. 
For given h > 0, let pi , p 2 be the respective minimizers of the following schemes: 



Pi := argmm 

pev 2 ,M z (R d ) 



Em\p\+yWi{p,PU 



fori = 1,2, 



(1.9) 



Then p\ < P2 a.e.. 
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This comparison result is new in the context of Wasserstein distances and might be of independent 



interest (see section 5.1 for more discussions). As a consequence one obtains geometric properties of the 



discrete solutions such as the confinement property (Corollary 5.5 1 



Lastly, making use of this confinement result, for strictly convex $ (but not necessarily uniformly 
convex), we have the following result concerning the long time behavior of p^ starting from general 
initial data: 



Theorem 1.4 (Convergence to the stationary solution, sec Theorem 5.6 1. Let 2 < m < oo. Let <J? be 
strictly convex and satisfy (A2) and (A3'). Assume the initial data po S P2(R d ) has compact support, 
and in addition satisfies ||po ||oo <l™ the case m = oo. For 2 < m < oo, let p m be given as the gradient 
flow for E m with initial data po, as defined in Theorem \4-l\ b). Then as t — > oo, p m (-,t) converges to 
the unique global minimizer ps of E m exponentially fast in 2- Wasserstein distance. 

1.3. An outline of the paper. In section 2 we introduce the notion of viscosity solutions for (P) and 



state basic properties of solutions. This part is largely parallel to IK] . In section 3 we show Theorem 3.5 
A key ingredient in this section is Theorem |3.4[ which identifies properties of the half-relaxed limits of 
p m as 77i — > oo. We point out that such convergence is previously known without the presence of the drift 
( [GQ| , [K]), but the presence of the drift and the resulting non-monotonicity of the support {p(-, t) > 0} 
causes new challenges. In particular the weak formulation used in |GQ] based on variational inequalities 
no longer applies, and thus we proceed with the viscosity solutions approach similar to those taken 
in [K]. The argument presented in Theorem |3.4| is of independent interest: it presents a strong stability 
argument which would apply to a general class of non-monotone free boundary problems. Let us point 
out that the assumption (Al) not only justifies (P) but also ensures the non-generacy of solutions of 



(P) near the free boundary which leads to stability properties (see e.g. the proof of Theorem 3.4 ) 



In section 4 we introduce the corresponding discrete-time schemes with free energy E m and E^ respec- 
tively, and we study the convergence of the discrete solutions (and continuous gradient flow solutions) 
as ?7i — > oo. There are new difficulties in handling the singular limit m — > oo, since the discrete solutions 



Pm.h corresponding to free energy (1.7 1 are not necessarily less than 1. Lemma 4.4 ensures that p m j t 
can be approximated with a density less than 1 which is close to the original solution in W% distance 
and has similar energy E m . This approximation as well as estimates between p m ,h and Poaji obtained 
in Proposition |4.5| enable us to prove Theorem |4.2| Finally, by combining the uniform convergence 
results obtained in Theorem 13.51 and Theorem 14.21 we conclude with Theorem 11.21 Let us mention that 



the T-convergence approach ( [DM| - Ser ) may apply here to derive the convergence of p m to poo in 
2- Wasserstein distance. On the other hand our approach is more quantitative and thus provides, for 
example, convergence rates in terms of m. 

Finally, in section 5, for any fixed 2 < m < oo, we present a comparison principle between solutions 



Pm.h of the discrete-time scheme corresponding to free energy E m when $ is semi-convex (Theorem 5.1 1. 
As mentioned above this result is new for the (discrete) gradient flow solutions in the setting of Wasser- 
stein distances. As applications of the comparison principle, we discuss some confinement results and 



the long time behavior of p m for convex $ in section 5.2 (Theorem 5.6) 



1.4. Remarks on possible extensions. For simplicity of the presentation we did not consider the 
most general setting our approach could handle. Below we discuss several situations where our approach 
(partially) extends. 

1. Our approach would apply, with little modification, to the problem confined in a domain £ C M. d 
with Neumann boundary data. On the other hand our approach would not apply, at least in 
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its direct form, if one puts an exit (e.g. Dirichlet) condition on parts of <9E. The challenge 
is in showing the convergence of discrete solutions, due to the fact that the A-convexity of the 
associated energy no longer holds. On the other hand, the analysis in section 2 and 3 should 



still go through to yield that the solution of (1.8 1 converges to the solution of (P) in domain E, 
with corresponding boundary conditions. 

In the case that A$ is not necessarily positive, and for general initial data < po < 1, the 



results in sections 4 and 5 are still valid and one can conclude that the solutions p m of (1.8) 
uniformly converges to a limiting profile p^ in 2-Wasserstein distance. In this case, the jammed 
region {poo{-,t) = 1} no longer satisfies finite speed of propagation and may nucleate at times. 



Due to this reason further characterization of p^ beyond as a weak solution of ( 1.5 ) remains open 



As mentioned above, without (Al), the unique characterization of the continuum limit p,^ given 
as the limit of discrete-time solutions remains open. In this case we suspect that a rather unstable 
mushy region {0 < p^ < 1} would develop in the limit m — > oo, generating non-uniqueness of 



Poo. A similar difficulty arises in the analysis of PQV where singular limits of degenerate 
reaction diffusion equations are considered. 

2. ON THE CONTINUUM SOLUTIONS 

In section 2 and 3 we assume that $ satisfies (Al). As mentioned before, we do not expect classical 



solutions to exist globally either for (P) or (1.8). Hence to investigate qualitative behavior of solutions we 
begin by introducing the notion of weak solutions for (P) , in our case the viscosity solutions. This notion 
of solutions is particularly useful when we are interested in the stability properties of interface problems. 
Let us point out that solutions of (P) may be discontinuous due to the quasi-static nature of the evolution, 
and due to the singularity of the free boundary. Therefore in the definition of viscosity solutions we need 
to consider semi-continuous functions, in contrast to the definitions of viscosity solutions in section 3. 
We introduce a definition using comparison with smooth functions similar to the one in IK] and IP]. 

Definition 2.1. A nonnegative upper-semicontinuous function u defined in Q := R d x [0,oo) is a 
viscosity subsolution of (P) with compactly supported initial data uq if the following hold: 



(a) u = u att — and {u > 0} = {u(x, t) > 0} n {t = 0}; 

(b) {u > 0} n {t < t} C {u > 0} n {t < t} for every r > ; 

(c) For every cf> 6 C 2 ' l (Q) that has a local maximum of u — <f> in {u > 0} D {t < to} at (xo,to), 

(i) if{x (h t )e{u>0}, -A0(a; o ,i o )< A$(x ). 
(ii) if (xq, to) G d{u > 0}, u(xo, to) — 0, and if \S7(fr(xo, io)| 7^ 0, then 

mm(-A<f)-A$,<t>t- |V</)| 2 - V0 • V$)(*o,*o) < 0. 

Note that the condition (c) (ii) is to ensure that limits of viscosity solutions are viscosity solutions, 
since the boundary can collapse in a limit and then boundary points of the limiting functions becomes 
interior points of the limit. 

Definition 2.2. A nonnegative lower-semicontinuous function v defined in Q :— R d x [0, oo) is a viscosity 
supersolution of (P) with initial data Vq if the following hold: 

(a) v = vq at t = 0. 
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(b) For every </> £ C 2 ' X (Q) that has a local minimum zero of v — <f> in M. d x (0,£o] a t (^ch^o)? 

(i) if(x ,t ) e{v> 0}, -Acf>(x ,t ) > A*(x ). 

(ii) if (x ,t ) G <9{v > 0},w(a;o,io) = and if 

|V^»(xo,*o)| 7^0 and {<f> > 0} n {v > 0}nB{x Q ,t ) ^ ® for some ball B (2.1) 

then 

max(-A0 - A$, <f> t - |V0| 2 - V0 • V$)(x , £ ) > 0. 



The condition (2.1) is to ensure that <f> touches v from below in a non-degenerate way. 



Let us define, for a function h in Q, the upper and lower semi-continuous envelopes of h: 

h*(x,t) := lim sup h(y,s), h*(x,t) := lira inf h(y,s). (2-2) 

e_>0 {|x-y|,|t-s|<e} e— S-0 { |a;— y| , 1 1— s| <e> 

Definition 2.3. u is a viscosity solution of (P) wit/i initial data uo z/u* and m* are respectively viscosity 
sub- and super solutions of (P) wii/i initial data uq. 

We will discuss several properties of viscosity solutions which will be used in the main theorem of the 
article. 

2.1. Inf- and Sup-convolutions. Next we introduce regularizations for viscosity solutions of (P), 
which is by now standard for free boundary problems (see e.g. [C V] ) . Given a viscosity subsolution u 
and r > 0, we define 

u r = sup u(y, t) for t > r (2-3) 

B r {x,t) 

and likewise given a viscosity supersolution v, and r, 6 > with (5 << r, we define 

v r = inf v(y,r) for £ > r. (2-4) 

B r _ 54 (x,t) 

These are called the sup- and inf- convolutions, respectively, and serve to smooth out viscosity solutions to 
help analyze the speed of the free boundary. The following properties of u r and v r are direct consequences 
of their definitions. 

Lemma 2.4. (a) u r is a viscosity subsolution of (P). Moreover, at each point (xo,£o) € d{u r > 0} 

there exists a space-time interior ball B such that 

B C {u r > 0} and~BC\{u r = 0} = {(x ,£o)}. 

(b) v r is a viscosity supersolution of (P). Moreover, at each point (xo,£o) G <9{w,, > 0} there exists 
a space-time exterior ball B such that 



B C {v r = 0} and B (1 {v r > 0} = {(x , to)}- 

Let e„+i denote the vector (0, ..., 1) in Q. The following two lemmas will prove useful in our analysis 
later. The first lemma can be proven with a parallel proof to that of Lemma 2.5 in IK] and thus we omit 
the proof. The second lemma is more interesting and involves ruling out the case of local total collapse 
of the solution, that is, the solution completely vanishing at a given time. The proof relics on (Al) to 
build a quadratic barrier subsolution. 

Lemma 2.5 {{u r > 0} cannot expand with infinite speed). Suppose (xo,£o) G d{u r > 0}. Then the 
corresponding interior ball cannot have its outward normal as e„ +1 at (:Eo,£o). 
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Lemma 2.6 ({v r > 0} cannot shrink with infinite speed). Suppose (xo,to) € d{v r > 0}. Then the 
corresponding exterior ball cannot have its outward normal as — e n +\ at (xo,to). 

Proof. 1. Suppose that {v r > 0} has an exterior ball with outward normal — e„ +1 at a point (xo,to). 
Then at (xo,ti), v will have an interior ball B\ centered at (xo,to) where t\ — to = r — 5to and B\ has 
outward normal e n +i at (xo,ti). 

2. Fix a number A satisfying 

A<- r^rr, A < 1, \«r-5t . 

5max B2(x„,i ) |V$| 

3. We define 

w(x, t) := v (Ax + x , A 2 (i - 1) + h) 
This serves to map the cylinder 

C := {(x,t) : |x-Xi| < A,ti - A 2 < t < t{\ 
to the cylinder C := {|x| < 1} x [0, 1]. Then to is a viscosity solution of a re-scaled version of (P): 

Aw + A 2 A$i =0 in {u > 0}; 

y = -d v uj - Aa„$i in d{u > 0}, 

where $1 is a rescaled and recentered version of $. By our choice of A, the bottom of C is strictly 
contained in B\, and so by lower semi-continuity we can find e > satisfying uj > e at t = 0. 

4. We construct our barrier. Define 

ip :=a(l-t/5-x 2 /2) 

where we choose a > so that a < min(infc A$i, e). Then — A<p — a < A$i, and on the bottom of 
C, (p < e < uj. On the sides of C, ip < < to, so ip < to on the parabolic boundary of C . However, 
ip(0, 1) = 4a/5 > = w(0, 1), so they eventually cross. 

5. We examine their crossing. To this end, we define T to be the first crossing time of (p and u>: 

T := inf{2 > 0| there exists x e C s.t. co(x,t) - <p(x,t) < 0} 
Then we can find a sequence (x n , t n ) with t n \.T and 

co{x n ,t n ) - ip(x n ,t n ) < 

Now we are in a compact set so we can suppose that i„->x6 C, and since ui is lower semi-continuous, 
we must have that 

u{x,T) - <p(x,T) = -p < 
Then (x,T) must be in the parabolic interior of C. 

The fact that this is a local minimum of u — (f> follows since it is the first time ui and <j) cross. We are 
done now because —Aip(x,T) — a < A$i(S,T) and 

(p t - \V(p\ 2 - Wtp ■ V$i = -a/5 - a 2 x 2 - \ax ■ V$i < 

where the final inequality comes from our assumption on A. □ 
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2.2. Comparison principle. The central property of the viscosity solution theory is in the comparison 
principle, which we state below. The proof is mostly parallel to that of IK] , and thus we only sketch the 
outline of the proof. 

We say two functions u, v : R — > R + are strictly separated, denoted by u -< v, if 



u < v in {u > 0} and {u > 0} is a compact subset of {v > 0}. 

Theorem 2.7. Let u and v be respectively viscosity sub- and supers olutions of (P). Ifu(-,0) -< v(-,0) 
then u(-, t) < v(-, t) for all t > 0. 

Sketch of the proof 

1. Due to the fact that uq -< vq, applying Definition 2.1 (a)-(b) and the semi-continuities of u and v, 
we have u r (-, r) -< v r (-, r) for sufficiently small r > 0. 

2. We claim that u r < v r for all times bigger than r, which yields our theorem. Hence suppose not, 
and define 

to := sup{t : u r (-, s) -< v r (-, s) for s < t} < oo. 

One can then proceed as in [K|, using the above lemmas to exclude the possibility that u r and v r cross 
over each other discontinuously in time, to show that at t = t , there is a point Xq such that 

(x , t ) € d{u r > 0} n d{v r > 0}. 

Moreover, there exists an interior ball Bi to {u r > 0} and an exterior ball B 2 to {v r > 0} at (xo,io) 
such that 

Let (i^TOi) be the interior normal to the interior ball B\ and (y^m-z) be the exterior normal to the 



exterior ball B 2 at (xo,to), with \v\ = 1. Due to the Lemmas |2.5| and 2.6 ?ni and rri2 are both finite. In 
particular at t = to both the sets {u r > 0} and {y_ r > 0} have the interior space ball B x fl {t = i } w hh 
interior normal ^. Since u r crosses t^. from below at (xo,to)j we have mi > m2- Moreover at (xo,io)> 
the support of u r propagates faster than normal velocity mi, and the support of v r slower than normal 
velocity mi- Formally speaking, we would like to claim from the fact that u r and v r are respectively 
sub- and (strict) supersolutions of (P) that 

|Vu r | > TOi + /i and |Vu r | < rn2 + (J,, where /x = V$(a;o, to) ■ v - (2-5) 

From the claim, we deduce a contradiction since u r < y_ r at t = to and m\ > mi. 



3. To prove (2.5 1 in the viscosity sense, we can use appropriate barriers to compare with u r and v r , 
to measure the growth of these functions at xq. This part of the proof is parallel to that of Theorem 
2.2 in K . Indeed the barriers corresponding to our problem (P) are constant multiples of the ones 
constructed in Appendix A of |KJ. □ 

Remark 2.8. Let us point out that, due to the restriction on the strict separation of the initial data, the 
above comparison principle does not immediately yield the uniqueness of the solutions for (P) . Later in 



the paper we will derive the uniqueness result (see Theorem 3.5 1, by showing that i^-contraction holds 



between the characteristic functions of the positive sets of the viscosity solutions. 
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3. Approximation by degenerate diffusion with drift 



As in section 2 we continue to assume (Al). Let p be a weak, continuous solution of (1.8), as given 
in [VJ. We define the pressure variable u by 

(3.1) 



m — 1 
Then u formally solves 

(PME-D) m u t = (m- l)u(Au + A$) + |Vu| 2 + V« • V$. 



In |CV| (for $ = 0) and in KL it was shown that u is a viscosity solution of (PME-D) m . For 



completeness we review the definitions. First we define a classical solution of (PME-D) m as a nonnegative 
function u € C 2,1 ({u > 0}) that 

(a) solves (PME-D) m in {u > 0}, 

(b) has a free boundary T — d{u > 0} which is a C 21 hypersurface, and 

(c) r evolves with the outer normal velocity |Vu| + n ■ V$, where rj is the inward normal of T. 

We then use the classical solutions as test functions to define viscosity solutions of (PME-D) m . 

Definition 3.1. A non-negative continuous function u defined in Q := R d x (0,oo) is a viscosity sub- 
solution of (PME-D) m if for every 4> € C 2,1 (Q) that has a local maximum zero of u — 4> in {t < to} at 

{<t> t - (m - 1)0(A0 + A*) - |V0| 2 - V0 • V*)(as 0) to) < 0. 
Definition 3.2. A continuous function v : Q — > R + is a viscosity supersolution of (PME-D) m if: 

(a) For every <p £ C 2,1 (Q) that has a local minimum zero of v — <p in {v > 0} D {t < t } at (xo,to), 

(4>t - (m - 1)0(A0 + A*) - | V0| 2 - V0 • V*) (a?o,*o) > 0. 

(b) Any classical solution of (PME-D) m that lies below v at time t\ > cannot cross v at a later 
time. 

Finally, u is a viscosity solution of (PME-D) m with compactly supported initial data Uq if it is both 
a viscosity subsolution and supersolution of (PME-D) m and both u(-,t) and {u(-,t) > 0} uniformly 
converge to uq and {uq > 0} as t — ¥ 0, respectively, in uniform norm and in Hausdorff distance. 

Let us point out that the above definitions, based on comparison with classical solutions, are essentially 
in the same spirit as the definition of viscosity solutions of (P) introduced in section 2. 

3.1. Properties of u m at the free boundary. We remark that the definition of viscosity supersolu- 
tions of (PME-D) m only applies in {v > 0} in order to make the viscosity solution notion be equivalent 
to the idea of weak solutions. This has the consequence of needing extra effort to analyze the behavior 
at the free boundary, which is provided by the following lemma. Its proof is analogous to Lemma 1.7 
in IK], with the difference in the construction of barriers. 

Lemma 3.3. Let v be a viscosity supersolution (subsolution) of (PME-D) m , and suppose that <f> is a 
smooth function where v — <fi has a local minimum (maximum) zero in {v > 0} at (xQ,t Q ) € d{v > 0} 



with to > 0. If 4> satisfies (2.1) at (xo,to), then 

(^-|V0| 2 -V0-V$)(x o ,t o )>(<) 0. (3.2) 
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Proof. First note that the subsolution case by definition is trivial as discussed above, since (f>(xo,to) = 0. 
Thus we proceed to the supersolution case. We may set to = after a translation. 

Let us fix constants r, 5 > and prove the lemma for the inf-convolution of v, 

W(x,t) = inf v(y,r) 

B r ._ st (x,t) 

Then the lemma follows by taking 6 — > and then r — >• 0. 



Now suppose that for a smooth 4>, W — <fr has a local minimum in {W > 0} at (xo, 0) € d{W > 0}, 
with <f) satisfying (2.1). By perturbing <f> we may assume that the minimum is strict. Let H be the 
hyperplane tangent to {<fi > 0} at (xo,0), with (z/, 7) the inward normal to H with \v\ = 1. Note that 
7 > —00 from Corollary 2.16 in [KLj . Let a — |V</)|(x o ,0) = 0„(a;o,O) > 0. Towards a contradiction, we 
assume that (3.2 1 fails, and so it follows that for some a > 

7 = Vf, = ^Oo, 0) < (o - <r) + 1/ • V*(x ). (3-3) 

Hence 7 is finite. Moreover we have 

W(x,t) > (j){x,t) in B v (x ) x [-77,0] for r\ « 1. (3.4) 

Due to the regularity of </>, there exists a space ball B Q interior to the set {x : 4>(x, 0) > 0} with 
x € <9i?o- We define 71 as follows: 

7 + cr/4 if7>V$(x )-^ 

7i ~ S 

a/2 + V$(a;o) • v otherwise 



Then we use the result of Lemma A.l to find a classical subsolution S of (PME—D) m in a neighborhood 
B ti (xq) x [—77,77] that firstly has initial support inside Bq, secondly has advancing speed 71 at (xo,0), 
and lastly has a parameter < e < min(cr/4, a/4) such that S satisfies 

7i > IVSI + z/ • V$ - e at (« ,0). (3.5) 

This condition helps us to show that it initially lies under <f>. 

We now claim that S lies under W in B v (xo) x [—77,0] for sufficiently small r\ , which will yield the 
desired contradiction to the fact that S is a subsolution and W is a supersolution, since S will cross W 
at (xq,0). 



Due to (3.3) and (3.5), we have 



if 7 > V$(a;o) • v 



\VS\(x ,0)<{ z (3.6) 

3a/4 otherwise 

<a=|V0|(zo,O). (3.7) 

On the other hand, observe that the support of S propagates with the normal speed faster than that 
of <p at (xo,0) due to (3.5). Due to the regularity of </> and S and their ordering at t = it then follows 
that 



{S > 0} C {(j) > 0} in B v (x ) x [-77, 0] (3.8) 

if 77 is sufficiently small. From the above two inequalities it follows that S < 4> in B v (xo) x [—77, 0] if r\ is 
sufficiently small. We can now conclude using the fact that (j) < W in that neighborhood. 
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□ 

3.2. Characterization of the half relaxed limits of u m as m — > oo. Let f^o and u$ as given in the 
introduction, and let u m be the unique viscosity solution to (PME-D) m with the initial data Uq. Recall 



that u m is given as the pressure variable of p m by (3.1), where p m assumes the corresponding initial 
data (^ao) 1 /'" 1 - 1 '. Let us then define 

ui{x,t) = inf sup u m (y,s); 

n>0 m >„ 

\(x,t)-(y,s)\<l/n 

u 2 (x,t) = sup inf u m (y,s). 



n>0 



rn>n 



\(x,t)-(y,a)\<l/n 

Note that the {u m } are uniformly bounded in m, as a consequence with comparison with stationary 
solutions of the form (C — $(x))+ with sufficiently large C > 0. Hence u\ and u 2 are both finite. 

Since we cannot guarantee that the support of u\ traces those of u m , we need to define an auxilliary 
function. Let us define the function 

i](x,t) := limsup X{supp(« m )}(y, a), 

m— >oo 
(j/,s)->(x,t) 

and the closure of the support of rj: 



ci = {0, s) -. r)(x, s) > o}, ci(t) -.= nn {s = t}. 

Finally, let us define the largest subsolution of the Poisson equation —Aw = A$ supported in fi: 
ui ■= [sup{w : R d x (0, oo) ->• R such that - Av < A$ and v = outside of CI}]*. 



Here /* denotes the upper semicontinuous envelope of /, as defined in (2.2 1. 

Note that then u\ = outside of CI and for each t > 0, U\ satisfies 

-ASi(-,*) < A$ in R d , -Aui(-,t) = A$ in the interior of Cl(t). 

This auxiliary function iii is indeed the new component of the proof compared to the corresponding 
theorem in [K]. We point out that u\ is positive in the interior of Cl(t) due to (Al). 

Theorem 3.4. Let u\, u 2 and u\ be as given above, and let $ satisfy (Al). Then ui is a viscosity 
subsolution of (P) in Q, and u 2 is a viscosity supersolution of (P) in Q with initial data uq. 

Proof. First note that u 2 is lower semicontinuous by its definition. Likewise, u\ is upper semicontinuous. 

A. u 2 is a supersolution: 

1. Suppose we have a smooth function <p and u 2 — (f> has a local minimum at (xq, to) in 
{u 2 > 0} n {t < to}. By adding e(t — to) — e{x — xq) 2 + c to <fi one may assume that the minimum is zero, 
and is strict in C r n {u 2 > 0}, where C r := B r (xo) x [to — r, to] for small r > 0. 

If (xo,to) is in {u 2 > 0}, by lower-semicontinuity of u 2 , we can make r smaller and assume that 



C r C {u 2 > 0}. On the other hand if (x ,t Q ) € d{u 2 > 0}, we can assume that (2.1 1 holds for 4>. In 
particular, |V0| 7^ so that U2 — > in C r n {w2 > 0} c away from (xq, to). Thus in either case we can 
find that u 2 — (j) has a strict local mininum zero in all of C r . 
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2. We now claim the following: if r is sufficiently small, along a subsequence u m — 4> has a minimum 
at points (x m ,t m ) G C r with (x m ,t m ) -> (#o,*o) an d {x m ,t m ) G {u m > 0}. 

To show this, define (x m , t m ) = argmin c (u m — </>); we can assume that the sequence only ranges over 
to that achieve the infimum of u 2 at (xo,to). Let (x 1 ,t') be a limit point of {(a; m ,i m )} m . 



First let us show that upon further refinement of our sequence we have (x m , t m ) G {u m > 0}. Clearly 
this is true if (xo,io) <= { u 2 > 0}, and thus suppose (xo,£o) S <9{u2 > 0} and (x m ,t m ) lies outside of 
the support of u m . Then we can assume that (2.1 1 holds for <f>, so in particular we can assume that 
|V</>| ^ in C r . This rules out the possibility that (x m ,t m ) lies in the interior of C r . Also, in this case 
(f>(x , to) — 0, and so we can find a > so that <f> < —a < on dC r n {u 2 > 0} c . This rules out the 
possibility that (x m ,t m ) lies on the boundary of C r . Thus we conclude that (x m ,t m ) £ {u m > 0} for 
sufficiently large m. 

Next let us verify that (x', t') = (xq, to). By definition for arbitrary (y, s) in C r 

(u m - 4>){y,s) > (u m - cf>)(x m ,t m ). (3.9) 

Since (x m ,t m ) — > (x',1/), for each n there is M(n) so that \{x mi t m ) — (x' ,t')\ < 1/n if to > M, and we 
may assume that M(n) > n. Then 

inf u m (x m ,t m ) > inf u m (y,s)> inf u m (y,s). 

m>M{n) m>M(n) m>n 

|(x',t')-(v,«)|<l/n |(x',t')-(j/,s)|<l/ri 

Taking sup ra on both sides we find limmf m _>oo u m (x m ,t m ) > U2{x',t'). Then, taking liminf of both 
sides of (3.9 1 as (y,s) — > (xo,to) and m — > oo, we find 

(u 2 - <f>)(x ,t ) > (u 2 - <f>)(x',t') 

which contradicts that (xo,to) is the strict minimum of u 2 — <t> in C r . This proves our claim. 

3. To finish showing that u 2 is a viscosity supersolution, take <j) and (x m ,t m ) as given above. When 
(xo,to) £ {u 2 > 0}, a straightforward computation using the the properties of u m as viscosity solutions 
of (PME-D) m gives 

-A(j>{x Q ,t Q )> A$(x ,t ) 

as needed. Next suppose (xq, to) £ d{u 2 > 0}, and that (2.1 1 holds for cf>. Suppose towards a contradic- 
tion that there is a > so that 

max(-A</> - A$, cf> t - |V0| 2 - V<f> ■ V$)(x , to) = -a < 0. (3.10) 



Let us define (f> m :— (f> + C(m) so that (u m — (f> m )(x m ,t m ) = 0. Since (x m , t m ) — > (xq, to), ( |3.10[ ) yields 
that 

{{<t> m )t ~ (m - l)0 m (A0 m + A*) - |V0 m | 2 - V0 m • V$)(a; m ,i m ) < 0, 
which contradicts with the fact that u m is viscosity solution of (PME-D) m . Thus we have (x m ,t m ) G 
d{u m > 0}. But then the inequality 

{{(f>m)t ~ |V0 m | 2 - V(t>m ■ V$)(x m ,t m ) < -a/2 < 0. 



contradicts Lemma 3.3 which applies since </> is smooth and so satisfies (2.1) at (x mi t m ) for large m. 

B. u\ is a subsolution 

The subsolution part of our theorem is harder to prove, since a smooth test function touching u\ at a 
free boundary point (a;o,to) from above in 17 cannot be extended smoothly to outside of fl so that the 
order is preserved. Thus the proof of B requires a careful study of the behavior of the free boundary 
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of Mi, which is achieved by studying the properties of Ui and u m . First note that Definition 2.1 (b) is 



satisfied due to Theorem B.l We proceed to show the property given in Definition 2.1 (c). 

1. It is straightforward from the definition of Cl(t) that 

{t*i(-,t)>0}Cfi(t). (3.11) 

Parallel arguments to the supersolution case yield that — Aiti < A$ in {ui > 0} in the viscosity sense. 
Thus it follows that 

ui < ui. (3.12) 

Suppose that we have a smooth function <f> and ui—<p has a strict local maximum at (xq, to) in f2n{£ < to}- 
As mentioned before u\ satisfies — Aui(-,t) < A<I>, indeed with equality in the interior of f2(i). Thus to 
check that u\ is a subsolution, it is enough to consider the case when Xq € dd(to) and Ui(xo,to) = 0. 



Note that in this case (3.12) yields that U\ — 4> also has a local maximum at (xo,to) in fl n {t < t }. 



Now suppose towards a contradiction that 

a := min(-A0 - A$, (f> t - | V</>| 2 - V<f> ■ V$)(ar , t ) > 0. 
Then since <\> and $ are smooth, it follows that for a small r > 

min(-A0 - A$, & - |Vc/>| 2 - V</> • V$) > 2a/3 in C r := B r {x ) x [t - r, t ]. 

2. Let r be the parabolic boundary of C r . We claim that T Oil C {(f) > Sq} for some ^o > 0. 

To see this, suppose the claim is false: this means that we can find (y, s) € T n Q fl {</> > 0} c . But then 
<Kz/i s ) — ano - so (^i — ^)(j/> s ) — ~4>{y^ s ) > which violates the assumption that u\ — 4> is strictly 
negative in T n C r . 

3. Now we proceed to show that u m < <f> on the relevant part of the parabolic boundary, that is, there 
exists some e > independent of m such that 

u m < 4> — e on r n supp(u m ) for sufficiently large m. (3.13) 



To show (3.13), suppose not. Then we can find (#&,£&) € rnsupp(w mfc ) where u mk {xk,tk) > (f>(xk,tk) — x, 
and by compactness we can assume (xt,tk) — > (x',1/) € T D f2. Then we have that for each n, there is 
K(n) so that |(a;fc,ifc) — (x',t')\ < 1/n and mj. > k if fc > K(n), where we can assume K(n) > n. Then 

sup u mk (x k ,tk) < sup u k (y,s)< sup u k (y,s). 

k>K(n) k>K{n) k>n 

\(x' ,t')-(y,s)\<l/n \(x',t')-{y,s)\<l/n 

Taking the infimum over both sides, we find 

Ux{x',t') > lira sup u mk (x k ,t k ) > lira sup 4>(x k ,t k ) = 4>{x\t') > S , 

which contradicts that Ui — < on T n fi. 

4. Now we define 

£(»,*) = 0(^ + 7^,t), where v = -— -y^o, t )- 

|V0| 

Here 7 > is chosen small enough to satisfy first that 

min(-A£ - A$, & - |V£| 2 - V£ • V$) > 2a/3 + 0( 7 ) > a/3 in C r . 

and secondly that for m large, u m < £ on T n supp(u m ) (which is possible since £— 4> = 0(7) and w m — 
is bounded away from zero on T D supp(u m )). 
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(^O)^o) 




t = t 







t = t -r 



d{£ > 0} 



x 



Figure 3. The motivation for £: it crosses ui at an earlier time 



This justifies the following definition: 



T m '■— sup{£ : (u m — ^)(x,t) < for all x € {u m (-,t) > 0} fl C r }. 

Then r m will be the first crossing time of u m and <f>, provided they cross (since u m is continuous, we 
need not worry about jumps inside its support). 

5. We now wish to show that, along a subsequence, u m crosses £ in C r . To do this, we first prove that 
there is a subsequence {m^} so that 

C s C\{t < t } n supp(w m J ^ for all 6 > 0. 



To show this first observe that, since Mi(a;o, to) = 0, there exists M so that 



sup 

m>M 
\(x ,y )-(y,s)\<l/AI 



u m (y,s) < 1. 



(3.14) 



Now assume towards a contradictio n tha t there is a S where our claim fails; we can take S small enough 
so that 8 < M~ 1 . We use Theorem B.l to derive our contradiction. First, we use the theorem to find 



positive numbers r max , T that depend on K = 1, the behavior of $ near (xo, to), and dimension, and we 
may assume T < 5/8. Now set ro a number smaller than min(r max , 5/8) and /i = min(T/8, ro/8). By 
definition of O, we can find (x',t') within distance fi of (xo,to) where r](x',t) = 1. Thus we can find a 
subsequence {mk}'^ =1 all bigger than M and points (y mfc , SmJ € supp u mh within distance /i of {x' , t'). 

Consider a specific m^. Then by the assumption that u rrlk = in Cs fl {t < to}, since tq and Tq are 



chosen much smaller than 5, we find u mi = in B ro (x') n {i = to — T/2}. Further, by (3.14) 



u m fc < 1 on the parabolic boundary of -B2r (a; ) x [t , t + T] 



Thus we apply Theorem B.l to find that 

Hn k = G B ro /4(x') x [t - T/2, t + T/2] 



II,, 
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and by the size of tq, we have that y mk € B rQ / 4 (x') and s mk < t + T/2. This yields u mk = in a 
neighborhood of (y mfc j s m k )- This is a contradiction, so we find our claim holds for every to&, giving us 
our desired subsequence. 

6. We will use the subsequence from the previous step to show that r mk < to. Indeed note that, for 

|a|,/3>0, 

£(x + a, t - 0) = -7|V(/>(x ,fo)| +«■ V0(x o ,io) 

-/3&(a:o, to) + 0(|a| 2 + P 2 + a/3 + 7 2 ). 

Thus there exists S = 5(j) > so that if \a\ 2 + (3 2 < S 2 then £(xo + a, to — j3) < 0. Due to the previous 
claim we can now find points (y mh , s m J € Ca n {t < to} fl supp(u mfc ), hence 

(«m fc -0(Vm h ,S mh ) > 0. 

Thus r mjt < S mk < t . 

7. Consequently there exists a crossing point {x mk ,T mk ) <E C r C\ {t < to} n supp(u mfc ) where (u mk — 
0( I m*)'mJ = 0. Further, since u mk - £ < on T n supp (u m J from step 3, we have that (x mfc ,r m J 
is on the parabolic interior of C r . Then we have that 

min(-A£ - A$,6 - |V£| 2 - V^ • V$)(x mfc ,r m J > a/3 

which forces that (x mkl T mk ) £ d{u mk > 0}. But then the inequality 

[6-|V£| 2 -V£-V$](x mfc ,T m J>0 



contradicts Lemma 3.3 which applies since £ is smooth and thus satisfies (2.1) at (x mk ,T, 



"mjfe) imi; 



C. ui , U2 converges to uo at t = 

It is not hard to check via comparison with radial barriers of (PME) m , based on the local Lipschitz 
geometry of dilo, that Q(t) and {u 2 (-,t) > 0} converges to f2 m Hausdorff distance as t — > + . From 
this fact and that u solves — Au = A<I> in the interior of Oq, we have lim T _i.o ^i('j t ) — u o from the 
definition of u\. On the other hand u 2 satisfies — Ait2 > A$ in {u 2 > 0} n {t > 0} and thus we have 
liminf T _j.o u 2(', t) > uq. Since u 2 < Wi by definition, it follows that u 2 {-, r) converges to uq as r — > as 
well. D 

3.3. Convergence of u m asm-) oo. Now let us fix a compact set £l in M. d with Lipschitz boundary, 



and let uq be as given in (1.1). Let u m be the viscosity solution of (PME) m with initial data Uq. If 



we knew that {u m } locally uniformly converges to a function u as m — > oo, then Theorem 3.4 would 
yield that u is a viscosity solution of (P). Unfortunately we do not know whether such convergence is 
true: due to the quasi-static nature of (P), u may not be continuous over time and this may complicate 
the convergence of u m . Thus we take the alternative approach to show the convergence of the support 



of {u m > 0} (see Theorem 3.5 (b)). The proof relies on the fact that {u m } has a stability property 



obtained from the L 1 contraction of the corresponding density function p m given by (3.1). Using this 



stability as well as the comparison principle (Theorem 2.7) we will obtain the support of u m converges 



to that of the unique solution u of (P) . From this result we then obtain the uniform convergence of p„ 



to the characteristic function of fl t away from the boundary of fit (Corollary 3.6). 

Theorem 3.5. Take Oo an d uq as given above and let $ satisfy {Al). Then the following hold: 

(a) There exists a unique evolution of compact sets {£lt}t>o such that any viscosity solution u of (P) 
satisfies fit — |u(-,t) > 0} for each t > 0. 
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(b) For each t > 0, the Hausdorff distance djj(Q tl {u m (-,t) > 0}) goes to zero as m —¥ oo, and 
limsvL'p m ^ roo u rn (-,t) is uniformly bounded. 

Proof. 1. The proof is based on the L 1 -contraction property (see e.g. section 3.5 of (Vl), which states 
that for two weak solutions pi,p2 of (1.8 1, the L 1 norm of their differences decreases in time. In terms 
of the pressure variable pi 



rn-l"i 



this reads 



\m 



l/(m-l) l/(m-l) 



P-2 



')(-.*)llLi(R-)<ll(pi 



l/(m-l) l/(m-l) 



Pi 



li-Mwy 



(3.15) 



Let us fix the initial data uq and vq so that uq -< vq. For each m, let u m and v m be respectively 
the viscosity solutions of (PME-D) m with their respective initial data uq and vq. Let us consider u\, v,2, 
ui as given in Theorem |3.4| and let Viji^j^i denote the corresponding functions given in Theorem |3.4 
defined with {v m } instead of {u m }. 



Since uq -< Vo, Theorem |2.7| applies to u\ and i>2, and so using Theorem 3.4 yields that 

U\ < U\ < V2- 



On the other hand, (3.151 yields that 

The above inequality and the fact that u m < v m and u\ < v?. imply that 
\{v 2 (;t) > 0}-{«i(-,t) > 0}| < limsup||(^( m - 1 )-^^ 1 ))(-,i)|| Ll < \{v > 0}-{u > 0}|. (3.16) 

2. Take uq as given above, and let us consider 

V(x,t) := (inf{t> : v is a viscosity supersolution of (P) with uo -< u(-,0)})* 

and 

U(x, t) :— sup{u : u is a viscosity subsolution of (P) with u(-, 0) -< u }. 

Here /* denotes the lower semicontinuous envelop of /, as defined in (2.2). Due to Theorem 2.7 U 
(V) then has the property of being below (above) any viscosity supersolution (subsolution) of (P) with 
initial data Uq. 



Let us consider a sequence of initial data v n and Vq n such that 

(a) Vg' n -< u Q -< Vq U for each n; 

(b) Uq n uniformly converges to u and {vq n > 0} converges {uq > 0} uniformly in Hausdorff 
distance. 



±,n 



±,7l 



Such Vq n can be constructed using the fact that d£lo is locally Lipschitz. Now let uf' n and u^' n be 
the corresponding versions of u\ and U2 with the initial data v ' . Then due to Theorem 
definition of U and V we have 



3.4 



and the 



~i 1 '" < [/, V < u 2 '" for any n. 



Using these approximations of initial data, the fact that {V(-,t) > 0} is open, and (3.161, we conclude 
that 

(3.17) 



fit :={V(;t)>0} = {U{; t)>0}. 
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Now for any viscosity solution u of (P) with initial data u 0l we have U < u < V. Thus {u(-, t) > 0} = f2 t , 
and we showed (a). 



3. By Theorem 3.4 and the definition of U and V, we have 

U <u 2 < (ui)* < V. 
Hence we have 



n t = {«2(-,t)>o} = {tti(-,t)>o}. 

The above inequality and the fact that u% is a viscosity subsolution of (P) with initial data uq yield 

(b). □ 

In terms of p m — (—Mm) 1 "" 1 " 1 ' the convergence results can be stated as follows: 



Corollary 3.6. Let (flt)t>o t> e the family of compact sets in K as given in Theorem 3.5 and let p, 



solve ||1.8|| with initial data p m (-,0) = ( 2Z ^7 i Uo) 1 ^ m_1 ^ Then for each t > 0, 



(a) limsuppm < 1; 

(b) {p m (-,t) > 0} uniformly converges to fl t in Hausdorff distance; 

(c) p m (',t) locally uniformly converges to 1 in Int(il t ), and to in (il t ) c . 

The same result holds for p m with initial data \Q ■ 
This concludes our analysis on the limiting profile of p m . In the next two sections we study the 



gradient flow solution p^ of the crowd transport equation (1.4 1. Among other things, we show that p m 
converge to poo as m — > oo in the Wasserstein distance (see Theorem 4.2), and hence p^ must coincide 
with \n t ■ 

4. Convergence of the gradient flow solution as m — > oo 

4.1. Definition of the gradient flow solution and the discrete scheme. For section 4 we introduce 
more assumptions: 

(A2) inf $ is finite, and without loss of generality we assume inf $ = 0. 

(A3) $ is semi-convex, i.e. there exists AeK such that D 2 Q(x) > X(Id)dxd f° r all X € R d . 

(A3') In addition to (A3), HA^H^ < C for some finite C. 

The semi-convexity assumption (A3) guarantees the well-posedness of the discrete-time JKO solution. 
When we prove convergence results as to — > oo, we will replace (A3) by the stronger assumption (A3'). 
It ensures that A$(a;) cannot be too large, which makes it possible for us to obtain some quantitative 
estimates on the difference between p m and p^ for large to. (A2) is a technical assumption, and will 
be used explicitly in the proof of Lemma 4.3 in section B] The assumption (Al) will be only used in 



section 4 to link p^ with the free boundary problem (P), and in section 5 to obtain convergence results 
as t — > oo. 

We denote by ^(R ) the space of Borel probability measures on K d with finite second moment, i.e. 
the set of probability measures p(x) such that J Rd p(x)\x\ 2 dx < oo. For a probability density p £ "P 2 (IR d ) 
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we define its "free energy" E m [p] as 



1!) 



E m [p] :— S m [p] + / p{x)<&{x)dx for 1 < m < oo, 



(4.1) 



where j Rd p(x)$(x)dx corresponds to the potential energy of p, and 5 m [p] is its "internal energy", given 

by 



5. m [p] := / — p m (x)dx for 1 < m < oo, 
Jrj m 



while 5*00 is defined as 



Soo[p] 



for HpII^oo^) < 1 
-oo otherwise. 



(4.2) 



(4.3) 



Next let us introd uce the following discrete-time scheme (also called a minimizing movement scheme) 

We consider a time-step h > 0, and initial data po G P20R. ) satisfying ||po||.L~(R<i) < 1- 



JKO 



introduced by 

For 1 < m < oo, the sequence (p^ k ) n6 ti is recursively defined by p® m h — Po and 



1 
Zh 



Pltl e arginf \ E m [p] + ~W*(p^ h ,p) : p G ^ 2 (M d ) 



(4.4) 



where W^(-, •) is the 2-Wasserstein distance (for definition, see e.g. AGS .) We then define p m .h{ x , t) as 
a function piecewise constant in time, given by 



Pm,h{ x > *) : = Pm fi( x ) for t e I n/l > ( n + l ) h ')- 



(4.5) 



Under the assumption in (A3) that $ is semi-convex, one can check that for all m > 1, the free energy 
E m [p] is A-convex along the generalized geodesies with respect to 2-Wasserstein distance, where A is as 
given in (A3) (For the definition of generalized geodesies and A-convexity, we refer to Appendix O). 
One can then apply the theory of gradient flow solution developed in AGS , which gives the following 



existence and uniqueness results of the discrete solution, as well as a convergence result as h — ¥ 0. 



Theorem 4.1 ( AGS] ). Let 1 < m < oo and suppose $ satisfies (A3). Moreover suppose E m [p ] < oo, 
where E m be as given in (4.1 ). Then for given h > the following holds for the sequence (p 7 ^ h ) n <£N cls 
defined in (4.4 1; 



(a) Existence & Uniqueness for discrete solutions (Section 2-3 of JAGSf ): Let A be as defined in (A3), 
and let ho = — t for A < 0, ho — oo for A > 0. Then for < h < ho, p 1 ^ h is uniquely defined for 
all n eN. 



(b) Uniform convergence as h — > (Theorem 4-0.7 - 4-0-10 in \AG9p : Assume that $ satisfies (A2) 
in addition to (A3), and consider initial data po such that E m [po] < M for some constant M. Let 
p m ,h be as defined in (4.5). Then for any T > and step size < h < 1, there exists some p m (t) 
(and Poo(t) in the case m = ooj in Cv^([0,T],'p2(K d )) such that 

W 2 {p m , h {-,t),p m {-,t)) < C(X)VMhe- XT for all t e [0,T], 

where A is given by (A3). Here we say p G C w ([0, T];V 2 (R d )) if p(-,t) eV 2 {R d ) for each < t <T 
and 

p(-,t) ->/»(-, t ) weakly in V 2 {R d ) as t ->■ t in [0,T]. 



Moreover, for finite m, p m (x,t) coincides with the viscosity solution of (1.8 1 
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(c) Contraction in Wasserstein distance (Theorem 4-0-4 (iv) in \AG3f ): For a given m, consider the 
initial data poi,Po2 £ "Psffi- )> with E m [poi] < oo for i = 1,2. Let pi(x,t) and p 2 (x,t) denote the 
limit solutions as defined in part (b), with initial data poi o,nd po2 respectively. Then we have the 
following stability result, where X is as given in (A3): 

W 2 {pi{;t),p 2 {;t)) < e- xt W 2 (poi,p Q 2) for all t > 0. 



The above theorem yields the gradient flow solutions p m (-,£) and Poo(-,t). In this section, our main 
goal is to prove that asm-) oo, p m (-,t) converges to Poc{-,t) uniformly in t e [0, T] in 2- Wasserstein 
distance. Convergence rates will also be obtained in terms of m. Although the rate is not optimal, to the 
best of our knowledge, our result is the first that gives some explicit convergence rate as the exponent 
m — > oo in the porous medium equation. More precisely, our main theorem in this section is as follows: 

Theorem 4.2. Let $ satisfy (A2) and (A3'), and consider p 6 V 2 (M. ) satisfying ||po||oo < 1 o,nd 
J po&dx < M. Let p m (t) and Poo(£) be as given in Theorem^- lVb) with the initial data po. Then for 
any T > 0, we have 

lim sup W 2 (p m {t), Poo (<) ) = 0. 
m ->-°°te[o,T] 

More precisely, we have the following convergence rate: 

sup W 2 {p m (t),Poo(t)) < 
te[0,T] 



C(M,T, ||A$|j 
mV24 



We point out that under the additional assumption (Al) and the assumption that po = Xsi j we can 
combine the results in Theorem 4.2 with Theorem 3.5 and immediately obtain that poo must coincide 

Without these two additional assumptions, 



1.2 



with the free boundary 



with xn t almost everywhere, which gives Theorem 

Theorem |4.2| still holds, but our approach fails to yield the connection between p Q 

problem (P). Thus a further characterization of Poo beyond as a weak solution of (1.5) remains open in 

the general context. 



The rest of this section will be devoted to proving Theorem 4.2 In section 4.2 we consider the discrete 
JKO scheme (4.4) for E m and E^ respectively, with the same initial data ||po ||oo < 1- We show that if 



we run the JKO scheme for one step only, then their Wasserstein distance is small. Once we have the 
one-step estimate, we are finally ready to prove Theorem |4.2| in section |4.3[ which says the Wasserstein 
distance between the continuous gradient flow solutions p m and p^ also goes to zero as m — > oo, with 
an explicit rate in terms of m. 



4.2. One-step estimate for large m. We consider the initial data po G 7^2 (K rf ) satisfying ||po||oo < 1 
and with finite potential energy, and let h be some fixed small time step. Then for any 2 < m < oo, we 
define p m (and p^ in the case m = oo) as follows: 



Mr. 



argmm 

pe-p 2 (R d ) 



E m [p] 



1 
2h 



wi{ Pa ,p) 



(4.6) 



Our main result in this subsection is Proposition |4.5[ which says that the Wasserstein distance between 
p m and poo is of order 0(to -1 / 8 ) for large m. To show that we first establish the following two technical 
lemmas concerning p m for 2 < m < oo. 

Lemma 4.3. Let 2 < m < oo, and let $ satisfy (A2) and (A3), and consider the initial data po 6 
V 2 (JL ) satisfying ||po||oo < 1 a nd f p$ < M. Letting p m be defined as in (4.6), the following estimate 
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holds (where a + := max{o, 0}): 



/ (Mm - l)+dx < 2\ 



M + l 



Proof. Our proof is based on the following crude estimate: ^- J Rd (p m ) m dx < M + 1. This inequality 
directly comes from the fact that E m [p m ) < E m [po], together with the assumption (A2) that inf $ > 0: 

-(p m ) m dx < f p <S>dx + f -{po) m dx <M + 1, (4.7) 

which upon rearranging gives 

/ ( f i m ) m dx<m(M +1). (4.8) 

JWL d 

Note that for to > 2, we have 



' (p m ) m dx> f (l + m( Pm -l)+ m( - m 1 \ l i m -l) 2 )dx. 



(4.9) 



Combining the inequalities (4.8) and (J4.9I) together, we have 



x9 , 2(M + 1) 4(M + 1) 

i 771 — 1 to 

Finally, note that |{a; : p m (x) > 1}| < 1, and so the Cauchy-Schwarz inequality yields that 



/ (Mm- 1)- 



dx<2\l — . D 



The following lemma says that for large to, we can find a probability density fj, m that is close to \x m in 
Wasserstein distance, has maximum density bounded by one, and has potential energy not much larger 
than p m . 



Lemma 4.4. Let $ satisfy (A2) and (A3'). Under the conditions of LemmaUjft there exists a proba- 
bility density ji m £ ^(R ), such that ||Mm||L°°(R d ) < 1> 



/ /} m $diz;< f /i m $ C fe + 2||A$|j 00 W^-^, (4.10) 

./»<* jR-i V TO 

and /j m is "close" to p m in the sense that 

2(M 4- I") 1 / 4 

W 2 (/w m )< l X/J ■ ( 4 - n ) 

TO, 1 / 4 



Proof. Due to the previous lemma, J Rd (p m — l)+dx < 2J M ^ 1 for all 2 < to < oo and /i > 0. We denote 



by a := 2\J ™+ 1 for short, and note that a is small for large to. Next we will give a explicit construction 
of p m , such that it satisfies all the requirements. 

We begin with breaking p m into the sum 

p m (x) = pl n {x)+ p 2 m {x), 

where 

p] n (x) := mm.{p m {x), 1 - a}, p 2 m (x) := (/i ra (i) - (1 - o)) + . 
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The idea is to construct \x m by keeping /i^ and modifying \J? m . We first make the observation that 



[i m only contains a small amount of mass: more precisely, 



p? m (x)dx < la. 



(4.12) 



This is due to the following two facts. First, due to Lemma |4.3[ the mass of \x m above 1 cannot 
exceed a. Second, we claim |{/i m > 1 — a}\ < 1. To show the claim, suppose not, then we have 
J min{^t m , 1 — a}dx > 1 — a. As a result, J fi m dx > (1 — a) + a = 1, where the (1 — a) corresponds to 
the mass below (1 — a), and a corresponds to the mass squeezed between (1 — a) and 1 due to our (false) 
assumption. 

Let us now construct jl m as follows: 

ji m (x):=n l m (x) + (g*n 2 m )(x), (4.13) 

where * denotes convolution and g(x) := ^Xb(o,r)i where R(d) is the dimensional constant chosen such 
that J Rd g(x)dx = 1. Note that although R(d) depends on d, we indeed have R(d) < 1 for all d > 1. 



We claim that fl m constructed in (4.131 satisfies all the requirements stated in the theorem. First 



note that the facts J g = 1 and g > imply that fl m is nonnegative and has the same mass as /z m . To 
show that ||/2 m ||oo < 1, it suffices to check \\g * ^mlloo < a - Since the mass of /tjti * s ^ ess than 2a, this 
inequality is a direct consequence of Young's inequality: 

\\g * Mmlloo < llMmlllllfflloo < 2a • - = a. 



Next we verify that the inequality (4.10) holds, which is equivalent to 

/ (9*l4n)$dx< / /J, 2 m $dx + HAfcHooO. 

This can be rewritten as 

ti 2 m {x) \{g * *)(») - *(»)! dx < IIA4IUC 



for all a; G 



Since /i^ has mass less than 2a, it suffices to show that 



I 



B(x,R{d)) 



^(y)dy - $(#) 



<^I|A*| 



where we used the fact that R(d) < 1 to get the right hand side, and note that UA'&Hoo is finite due 
to (A3'). The proof of this inequality is similar to the proof of the mean value property for harmonic 
functions, and hence is omitted here. 



Finally it remains to show (4.11), which is equivalent to 

W 2 {^ m ,P<m,) < 



2a. 



(4.14) 



We now heuristically describe a transport plan, which is not necessarily optimal. First, we keep the mass 
of ix] n at its original location, so that no transportation cost is induced. Second, for for every "particle" 
located at x in /i^, the transport plan is to distribute it evenly in the disk B(x,R(d)). (Again recall 
that R(d) < 1 for any dimension d > 1.) Since the mass of /x^ is no more than 2a, the total cost of the 
transportation plan is bounded by 2aR(d) 2 , which immediately implies (4.14). □ 



Now we are ready to state the following one-step estimation, which controls the Wasserstein distance 
between \i m and [i^: 



QUASI-STATIC EVOLUTION AND CONGESTED CROWD TRANSPORT 



23 



Proposition 4.5. Let $ satisfy (A2) and (A3'), and consider the initial data po € ^(M. ) with 
1 1 po 1 1 oo < 1 and J po$>dx < M. Let A be as given in (A3), and A~ := — min{0, A}. For any < h < 32 ^_ , 
let p m and p^ be as defined in (4.6| for the cases m finite and m = oo respectively. Then the following 
inequality holds: 

W 2 (p m ,p co )<^—C(M). 



Proof. Let us fix M and d Suppose the statement is false. Then for an arbitrarily large Aq > 0, there 
exist m > 2, < h < QnM ^n — such that 



32||A$|| 



W 2 (n m , Moo) = Am 1/s , where A > A Q . 

To get a contradiction, we will construct a new probability measure r\ G 7^2 (K^) with ||t?||j,o 
such that the following inequality holds if A is chosen to be sufficiently large: 



(4.15) 
*> < 1 



E m [v] + ^Wi( Po ,v) 



Eoo[v] + ^Wi(p ,v) 



< 



■^m [f^n 



i 

2ft 



W^(p ,Mm) 



£00 [^c 



1 

2li 



w2(p ,n o 

(4.16) 



This means that r\ would beat at least one of the minimizers in (4.6) for some m (m may either be finite 
or +oo), contradicting the definition of p m and poo. 

The prob ability density n is constructed as follows. Let p m be the probability density constructed in 
Lemma 4.4 and we denote by T m the optimal transport map such that (T m )#po — p m . Similarly, let 
Too be the optimal transport map such that (Too)#po — Poo- Then r\ is defined as 

»?=(^r TO + ir 00 )#po. (4.17) 

r\ is thus the midpoint between p m and p^ on their generalized geodesies, as defined in Sec 9.2 of [AGS] . 



Next we will prove that r\ satisfies the inequality (4.16). This is done by proving the inequalities 



<|4.18|)-(|4.20|): 

S m [rj\ +Soo[ri] < S,n [Pm } + Soo [Moo ] 



1 
m ' 



1 



n&dx < / p m $>dx 
1 



Poo^dx 



C(M)(||A*| 



A) 2A 2 A 



W 2 (po>»?) < 



1 
2h' 



III TO 1 / 4 

2 , 



Wi( Po ,p m ) + ttt-WUpq^oo) + 



C(M) - A 2 /8 



(4.18) 
(4.19) 
(4.20) 



If Aq is chosen to be a sufficiently large number depending only on M and d, since A > Aq, the sum of 
these three inequalities implies (4.16) (where we make use of the assumption that h < „^\- )■> thereby 
yielding a contradiction. 



To show (4.18), it suffices to prove that ||»j||c» < 1, since then S m [rf\ = — J n m dx < — and SooM = 0. 
Proposition 9.3.9 of |AGS| states that the LP norm with p > 1 is convex along the generalized geodesies, 
and thus 

2|M| P < ||£m||p + HjUoollp for all p > 1. 
Sending p —$■ 00 in the above inequality yields ||??||oc < L since both \\p n 
byl. 



and 



are bounded 



(4.19 1 comes from the semi-convexity of $ given by (A3) (which is a consequence of (A3')). Let AeR 
be as given in (A3). Proposition 9.3.2 in AGS yields that J p&dx is a A-convex functional of p along 
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any generalized geodesic, and thus 

2 r]<S>dx < p m $dx + \ Hoo$dx- -AWi(fi m ,IJ>oo) 

JR d JR d JwL d 4 

GWIIAilU 



< 






(4.21) 



Poo^dx 



,1/2 



-XW 2 {p m ,p OQ ), 



where the last line comes from (4.10). Next let us estimate W 2 (MmiMoo)- Due to our assumption (4.15) 



in the beginning of this proof and the inequality (4.11), we have 



W 2 {fim,Hoz) < W 2 (n m ,Hoo) + W 2 (p m ,p m ) < Am' 1 ' 8 + C{M)m- l ' A . 
We take the square of the above inequality and apply the inequality (a + b) 2 < 2a 2 + 2b 2 , and obtain 

WK/i^/ioo) < 2A 2 m- 1 ' A + CiMjm- 1 / 2 . 



Plugging this inequality into (4.21) yields (4.19) 



Finally it remains to show (4.20 ). Due to Lemma 9.2.1 of AGS W 2 (po, •) is 1-convex along generalized 
geodesies, and thus 



Wi( P0 ,v) < jWaO*). M + 2^2 (PO, Moo) - -W 2 2 (/im,Moo)- 

Now, by the triangle inequality, we have 

1 / \ 2 1 1 / 

W 2 (po,TJ) < -\W 2 (po,fi m ) +W 2 (Mm,Am)J + ~z W 2 (/9 > Moo ) - T (^ W 2 (/Urn , Moo ) ~ W 2 {p m ,p m ) 

< 2^l(P0,Mm) + 2^2 2 (P0,Mco) " ^Wf (Mm, Moo) + ^JJ (1 + VF 2 (^ m , Poo)) 

lrrr2/ n2 l T r,2/ s C(M) - A 2 /8 

< -^|(p ,Mm) 2 +2^2(P0,Moo)+ V ^ 1/4 ' ■ 



(4.22) 



4.4 



as well 



For the second inequality we used the fact that W 2 (p, m , p m ) < C(M )m -1 ' 4 due to Lemma 
as that W(po, p m ) < C(M) for all m and h (otherwise p m would fail to be a minimizer of *f> m )- For the 
third inequality we use the assumption that W^MmiMoo) = AmT 1 '^. Finally, dividing both sides of the 



above inequality by h yields (4.20) 



□ 



4.3. Convergence of the continuum solutions as m — > oo. In this subsection, we give a proof of 
Theorem |4.2| The proof is done by combining the one-step estimation results in section |4.2| with the 
convergence results for discrete solutions as h — > 0. 



Proof of Theorem \4.3\ 1. Note that the assumptions on p immediately imply that E m [p ] < M + 1 
for all 2 < m < oo. This enables us to apply Theorem |4.1[ b): For all time steps h satisfying < h < ho 
(where ho is a small constant depending on M,T, ||A$||oo), we have 

W 2 (pm(t),Pm,h(t)) < C{\)VMe- XT Vh =: CVh for all t e [0,T], (4.23) 

and this inequality holds for both finite m and m = oo. For notational simplicity, the various constants 
C appearing in this proof may depend on M, || A^Hoo, T, and the value of C may differ from line to line. 

2. Now we fix the small time step h such that < h < ho, and our goal is to show that 

W 2 (p m ,h(t),Poo,h(t)) < CVh for all t e [0,T] when m > C7i" 12 . (4.24) 
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Figure 4. Illustration of the tree structure. 



Before we prove this inequality, let us point out the proof is finished once we obtain this: by combining 



(4.24) with the inequality (4.23) (and note that (4.23) holds for both finite m case and m — oo), one 



immediately has 

W 2 {p m {t) 1 p 00 (t)) < cVh for all t € [0, T] given that m > Ch~ 12 , 
which concludes the proof and would give the rate W2(p m (t),p 00 (t)) < Cra~ x / 2i for t € [0, T]. 



3. To prove the inequality (4.24) in step 2, note that it is equivalent to prove 

W 2 (p" h ,p" fc ) < Cv^ for aU n < ^. 



(4.25) 



4.5 



From now on we will denote p 1 ^ h by p™ n (and denote p^ h by p 1 ^ ) for notational simplicity. Proposition 
then shows that W2(pl n , p^) is small for sufficiently large m. To deal with the case n > 1 we consider 



the tree structure as illustrated in Figure [4] 
Here for n > 2, r\ n is defined as below: 



argmin 

pev 2 (R d ) 



E m [p] + llh Wi{p^\p)\ for// >2. 



We point out that for n > 2, llpJ^Hoo < 1 holds by definition, and in addition we have j p\^ l <frdx < 
J po&dx < M. Hence by taking p 7 ^ 1 as the initial data, the one-step estimate in Proposition |4 . 5 1 yields 
that 

W 2 (pZ,,V n ) <S:= C(M)m- 1/8 for all n > 2. (4.26) 

Let us denote d n := W2{pm,p\] ), which satisfies 

d n < W 2 (pZ,V n ) + W 2 ( V n ,p^) < W 2 {p n mlV n ) + 6, (4.27) 
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where 8 is as defined in (4.26). Now it remains to control W^lPm'l")- Note that p 1 ^ and r\ n are 



minimizers given by the discrete-time scheme (4.4) with the same free energy functional E. m , but with 
different initial data p^ -1 and p^ 1 - To estimate d n in terms of d n -i, we use Lemma 4.2.4 of AGS 



E m [R n m ])] 



which states that the Wasserstein distance between two discrete solutions does not grow too fast. More 
precisely, it gives the following inequality 

Wi(p n m ,v") < e- 2X ~ h [WK/C-VST 1 ) + h(E m [p^] 

<e- 2X ' h {d 2 n _ 1 + ha n - l ), 

where A - := — min{A,0}, with A as given in (A3). We denote a n _i := E^p^ 1 ] 
satisfies the following properties: 



(4.28) 

^mfc], Which 



a„ > for all n e N+, and ^ a n < M + 1. 



(4.29) 



Finally, we plug (4.28) into (4.27) to obtain the following family of inequalities: 



rfi < 5 

d„ < 



-2\~h 



dl 



ha„ 



8 for n = 2, 3, 



(4.30) 



4. We next focus on the inequality (4.30), and our goal is to show that dr < C\fh for 8 sufficiently 



small (more precisely, 8 < ft. 3 / 2 would be enough). By taking the square of (4.30) and applying the 

h 



inequality 2ab < ha 2 + K-, we obtain 



d 2 n < (f + ft) e ^ 4A ^K 2 _ 1 + /ia„- 1 ) + (f + -)5 2 

h 

< (1 + Ch)d 2 n _ 1 + h{2a n - l +2h), 



(4.31) 



where in the last line we let 8 < h 3 ' 2 so that (1 + ^S 2 < 2ft, 2 . Also note that b n '■— 2a„ + 2ft, satisfies 



T,n=o b n < 2 ( M +T+1). Now by dividing by (1 + Ch) n on both sides of ( |4.31[ ) and summing the 
inequality from 2 to n, we obtain that 

n 

d 2 n <dj(l + ChY 1 - 1 + Y^ hb k (l + Ch) n - k for all n. 



fc=i 



T/h T/h, 



Hence as a result, we see that W-^/Om , p<x> ) = dr/h satisfies 

dn 



*T/h 



< 



,CT 



ft 3 + 2{M + T+ \)e CT h < CVh 



(4.32) 



as long as 8 < ft 3 ' 2 (recall that 8 = Cm 1 ' 8 , hence it is equivalent with m > h 12 ), and so we are 
done. □ 



5. Comparison principle and long-time behavior for gradient flow solutions 



5.1. Comparison principle for the discrete-time solutions. In the beginning of section |1J we 
have defined the discrete-time scheme (4.4) for the porous medium equation with drift (1.8). Since the 



comparison principle for the viscosity solutions of (1.8) is well-known (see e.g. |KL| ), it is natural to ask 
whether the comparison principle holds for the discrete-time solutions generated by (4.4) as well. In 



this section we prove that this is indeed true, but the proof is quite different from the continuous case. 
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In fact comparison principle-type results have been shown between discrete gradient flow solutions with 



L 2 distances, for instance in (CI, GK , etc. The novelty in our result is that we address the discrete 



gradient flow solutions with W% distances, for which nonlocal perturbation arguments are necessary. 

In order to define the scheme for two ordered initial data, we need to consider non-negative measures 
which do not necessarily integrate to 1. We denote by V2,A^ d ) the set of non-negative measures which 
integrate to A > and have finite second moment. We also generalize the Wasserstein distance Wi as 
follows: For two regular measures pi,p2 & F2,A(^ d ), we define Wi{p\,p-i) as 

Wi(pi,p 2 )~ inf f \T(x)-x\ 2 Pl (x)dx. 

Next we state the comparison result. 

Theorem 5.1. Let <I> satisfy (A3). For 2 < m < oo ; consider the two densities poi € T'2.M 1 (^ d ), 
P02 € r P2,M 2 (W 1 ) with the property M\ < Mq, and poi < P02 a.e. (In the case m — oo, we require in 
addition that ||/)oi||oo < 1 for i — 1,2). For given h > 0, let p\,p2 be the respective minimizers of the 
following schemes: 



argmin Fi{p) '■— argmin 

P&V 2 , M AR d ) pGP 2 .M,(K d ) 



E m [p] + ^W 2 (p,p Ql ) 



fori = 1,2. (5.1) 



Then p\ < P2 almost everywhere. 



Remark 5.2. The proof of above theorem does not directly use the semi-convexity of $, except to 
guarantee the existence of the pi . 

Before we prove Theorem |5.1| we first state and prove the following simple lemma, which can be 



informally stated as follows: Given that p\ is the minimizer for the discrete scheme in (5.1) and T\ is 
the optimal map between poi and p\ , if a part of /?oi is forced to be transferred by the map T\ , then T\ 
is still the optimal map for the rest of poi- 

Lemma 5.3. Let 2 < m < oo 7 and let h and poi be as given in Theorem\5.1\ We denote by p\ the 



minimizer of J-\ as given by (5.1), and let 7\ be the optimal mapping such that Tiffp m = p\. Consider 
an arbitrary function r\ : M d — > R such that < n{x) < 1 for all x £ R d , and let (p(x) := Ti#((l — r])poi). 
Then Tif£(i]poi) minimizes 

Hp)--= f (-^ + p) m + P^]dx + ^-W 2 {r 1 p m ,p) 
J Rd \m J 2h 

among all p£P 2 ^(M^), where M = J Rd r\p m dx. 



Proof. Suppose that the minimum of J 7 is achieved by another measure p £ V 2 m(K^), such that F(p) < 
■F(Ti#(?7/9 i)). We denote by T the optimal map such that T#(r/p) = p. The claim is then that we can 



find a better transfer plan of poi than p\ in (5.1), yielding a contradiction. We construct the transfer 
plan as follows. 
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First, we separate poi m to two parts: ?7Poi an d (1 — v)Poi- Then we use 7\ to push forward (1 — rj)p Q i, 
and use T to push forward rjpo\. The resulting measure would be equal to <p + p. Then it follows that 

Fi(<P + p)= [ (-{v+~p) m + {v+~p)®)dx+^-Wi{p 0ll { V + p)) 

{h# + p) m + & + p)$)dx + ~Wi(r)Poi,p) + ^-W? ((1 - v)(X>i,<P) 

lu/2m _„w ,^ ( 5 - 2 ) 

2/i 



< 



Hp)+ I V&das + ^Wiai-rfipot,?) 



< T{Ti#{vp)) + I <P*dx + ~Wf ((1 - r,)p 01 ,cp) 

j R d in 

= T 1 (T 1 #p)=T 1 (p 1 ) 1 



which contradicts the fact that pi is the minimizer of (5.1 1 and so we are done. □ 



Proof of Theorem \5.1\ First we point out that once we prove the comparison result for all 
2 < to < oo, it will be automatically true for the case to = oo as well, due to the one-step estimate in 
Proposition 4.5 To see this, let us denote by pi >m the minimizer pi when the free energy is E m . Then 
gives us pi^ m — ¥ pi,oo as to — > oo in Wasserstein distance. Therefore if we know that 



4.5 



Proposition 

pi,m < P2,m a.e. for all 2 < to < oo then it directly follows that /9i )00 < P2.00 a.e. 

Due to the above discussion, it suffices to prove the comparison principle for any fixed to satisfying 
2 < to < 00. Let Ti denote the optimal map such that Tj^poi = Pi for i = 1,2. We prove by contradiction 
and suppose f2 := {p\ > p 2 } has non-zero measure. We first claim that 

\T 1 -\n)\T 2 - 1 (n)\>o, (5.3) 

which directly follows from the inequality below: 

/ poidx= / pidx (since Ti#p i = Pi) 

> / pidx (from the definition of 51 and the assumption that |S1| > 0) 

J? (5.4) 

= / P02dx (since T 2 #p 2 = P2) 



I 



> / p idx (since p 01 < p 02 ). 



Let Cl s = {x e R d : pi(x) > p 2 (x) + <5}, and let A s = {x eR d : Pi{T-l(x)) < |,/92(T 2 (a;)) < §}. Since 
U,5 >0 J1 5 = fi and U S>0 A S = R d , ([O yields that 



(r i - 1 (fi 5 )nA 5 )\r 2 - 1 (fi) > for sufficiently small S > 0. (5.5) 

From now on we fix 6 such that the above inequality is true, and denote 

B:=(2r 1 (n«)ni4 a )\r 2 - 1 (n). 

By definition of the set B, it immediately follows that T\ maps B into the set where p 2 + S < p x < I, 
while T 2 maps B into the set where p\ < p 2 < i. (Note that these inequalities hold in the a.e. sense). 
These facts are illustrated in Figure [5] 
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FIGURE 5. Illustration of the set B,T\(B) and T 2 (B). Recall that Tt is the optimal 
map between poi and pi for i — 1,2. Moreover, the set B is chosen such that T\(B) C 
{\ > Pi > 92 + 6}, while T 2 (B) c { Pl < p 2 < £}■ 

Let p e :— C/OqiXb, where < e <C <5 is a sufficiently small number to be determined later. Let 



ifi := Ti$=(pqi — p c ) = (1 — e)/9iXTi(B)> then by applying Lemma 5.3 to the optimal plan T\ in comparison 
to T 2 , and subtracting — tp™ on both sides, we arrive at the following inequality: 



< 



-((pi + T^ Pc ) m - -p? ) dx + E[T X ] 
m m 

^(p 1 + T 2 #p e r--pT)<lr+E[T, . 
m m 



(5.6) 



where 



E[T t ] := J d ((T,#p £ )$ + ^m(x) - xfp^dx, i = 1,2. 



Next we state a simple algebraic inequality without proof. For all real numbers a and b satisfying 

(5.7) 



< b < a < j and m > 2, we have 



a m - l b < — (a + 6) m - — a m < a m - l b + Cb 2 , 

771 777 



where the constant C only depends on m and S. Using (5.7), (5.6) yields that 



f ^r 1 (7i#ft)<fe + £pi]< f V?- 1 (T2#Pe)dx + E[T 2 ]+ f C(m,6)e 2 p 2 2 dx 

JR d JR d JT 2 (B) 

< f p™-\T 2 #p t )dx + E[T 2 ] + C(m,5)e 2 (since p 2 < ] in T 2 (B)) 

JR d d 

< [ (%-\T*#p e )iix + E[T 2 ] + C(m, 8)e 2 (since Vl < Pl ) 

JR d 

< f pT- 1 (T 2 #p t )dx + E[T 2 ] + C(m,S)e 2 

jR d 



(5.8) 



where the last inequality holds since p\ < p 2 in T 2 (B), and supp(T 2 #p e ) C T 2 (B). 



.'(II 
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Similarly, we define ip 2 := T 2 #(pq 2 — p e ), and note that p 2 is the minimizer to (5.1). We then apply 



Lemma |5.3| to the optimal plan T 2 in comparison to T\ , and an argument parallel to that above yields 
the following inequality for ip 2 : 



I ^-\T 2 ^Pe)dx + E[T 2 ]< f V2 n - l (T 1 # Pe )dx + E[T 1 }+C(m,8)e 2 

JR d jR d 



< 



m — 1 
Pi 



(T^pjdx + E[T X ] + C(m, 5)e 2 (since ip 2 < p 2 ) 



(5.9) 



Note that in the set T 2 (B), p 2 is bounded above by i, hence ip 2 is just smaller than p 2 by order e in this 
set, namely p 2 < ip 2 + e/5 in T 2 (B). Combining this with the fact that the integral of T 2 #p t is also of 
order e, and the assumption that m > 2, we have the following: 



JR d JR d 



(|5.10|) provides us a link between the RHS of (|5.8|) and the LHS of (|5.9|), and so we arrive at 

2 



f v>T~HTi#Pe)dx< f p 2 n - l {T l #p e )dx + C{m,5y 

JR d JM d 



(5.10) 
(5.11) 



Next we show that (5.11 1 leads to a contradiction if e is chosen to be small enough. First, recall that 
t>i = (1 — e)piXTi(B): and p\ > pi + 8 in T\{B). Hence if we let e be sufficiently small, we would have 
k > P2 + | in Ti(B). Then we have 

\T 1 #p e )dx> 



,rn — l i 



{pi 



\m-l 



(Ti#p e )da; 



> 



> 



m — 1 i / v \ 

Pi +(g) 



771 — 1 



{Ti#p t )dx 
6, 



pr i m#p e )^+e(-)™- i ii P1 || L1(B) , 



which contradicts (5.11 1 when we hx <5 and let e be sufficiently small. This concludes the proof. 



(5.12) 



□ 



Remark 5.4. By sending the time step h — > 0, the comparison principle for discrete solutions imme- 
diately leads to a comparison principle for gradient flow solutions. Also, although we only prove the 
comparison principle for the energy J p m dx with 2 < m < oo, the proof can indeed be easily extended for 
1 < to < oo 7 and also the case when the entropy part is given by J p log pdx . 

5.2. Confinement result and long-time behavior. In this subsection, we show some applications of 
the comparison principle for discrete JKO solutions. The first application is the following confinement 
result for discrete solutions (hence continuous gradient flow solutions as well), given that $ — > +oo as 

\x\ — > oo. 

Corollary 5.5. Let 2 < to < oo and let <f>(x) satisfy (A3) and the additional assumption that 



lim 



\x\— »oo 



$(x) 



-oo. Assume the initial data po S £°°(R ) has compact support, and if m 



do 



we assume in addition that ||po||oo < 1. Then the support for the discrete solution p 7 ^ h will stay bounded 
for all n, where the bound of the support does not depend on n or h. 

Proof, o Case 1: 2 < m < oo. For any A > , let us lo ok for the global minimizer pa of the energy E. 



as defined by ^j\ among V 2 , A (R d )- Due to 

(=-i(C, -.(x» 



CJMTU Lemma 6], the global minimizer pa is given by 



PA 



1 
m—1 
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where Ca is chosen such that the total mass of pa is equal to A. Observe that for any A > 0, such pa 
is also a stationary solution for the discrete JKO scheme, and it has a compact support. 

Therefore for any po £ L°°(M. d ) with compact support, one can choose A to be sufficiently large such 



that po < Pa a.e. Then one can apply Theorem 5.1 and obtain that p™ h < pa &■&■, hence the support 



of pl l n h stays within the support of pa for all time steps. 

o Case 2: m = oo. In this case, we first point out that for any mass size A > 0, the global minimizer 
of Eoa among V 2 ,a is given by some characteristic function xs A , where Sa is the level set of the function 
$, i.e. 

S A = {x £ R d : *(») > C A }, 
and Ca is chosen so that xs A has mass A. Moreover, since xs A i s the global minimizer of E^ , it must 
be a stationary solution as well. 

Recall that po has compact support, ||/?o||oo < 1, and $(x) — > oo as |x| —> oo. Therefore if we let A be 
sufficiently large, we will have supp po C Sa, which implies that po < XS A - Since xs A is a stationary 
solution, the comparison result in Theorem 5.1 immediately implies that supp p^ h C Sa for all n and 
/i, and we are done. □ 

Lastly we briefly discuss the long time behavior of the gradient flow solution p m for 2 < m < oo, 
when <I> is strictly convex and bounded below in M. d . In this case, one can easily obtain that the global 
minimizer for E^ in Vi(K d ) is ps := XOi where O = {i £ ff 1 : $(ar) < C}, and C is chosen such that 
Xo has mass 1. 

Theorem 5.6. Let 2 < m < oo. Let $ &e strictly convex and satisfy (A2) and (A3'). Assume the 
initial data po £ ^(^ ) ftas compact support, and in addition satisfies ||po||oo < 1 OT i/ie case ttt. = oo. 
for 2 < m < oo, /ef p m &e given as the gradient flow for E m with initial data po, as defined in Theorem 
\4- l\ b)- Then as t — > oo, p m (-,t) converges to the unique global minimizer ps of E m exponentially fast 
in 2-Wasserstein distance. 

Proof. If $ is uniformly convex in M. d , then there exists some A > 0, such that D 2 <fr(x) > XI for all 



x £ M. . In this case we can directly apply the contraction result in Theorem 4.1 c) between p 00 {x 1 t) 



and ps(x) (where ps is the global minimizer for the free energy E m [p] in T'2(^ d )), which gives 

W 2 (p oo (;t),p s (-))<W 2 (p ,Ps)e- Xt , 

and hence the 2-Wasserstein distance between p oc (a;,i) and ps{x) decays exponentially fast in t. 

On the other hand, if $ is strictly convex in M. d but not uniformly convex, we will make use of the 



confinement result in Corollary 5.5 As long as po is compactly supported, the proof of Corollary |5.5 
shows the support of Poo(-,t) will stay in some compact set Oa for all time, and indeed one can find an 
O a such that it is independent of m for all 2 < m < oo. This confinement result allows us to apply the 



contraction result in Theorem 4.1 (c), which gives that 

W 2 ( Poo (;t),p s (-)) < W 2 (p ,Ps)e-~ Xt , 
where A = inf{A : D 2 Q(x) > XI for all x £ Oa] is a strictly positive constant depending on po and 



Finally we remark that for finite m and p m , the corresponding result is shown in CJMTU , where 
they use entropy dissipation methods. We suspect the convergence rate to be exponential in stronger 
norms instead of Wasserstein distance, but this issue is not pursued here. 
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Appendix A. Constructing a (PME-D) m subsolution 

Lemma A.l. Fix e > 0, m > 0, a number -f, a point x' , and vector ft. Then if 7 > V$(V) • ft, 
there exists a positive constant r\ depending on e so that we can construct a classical subsolution S of 
(PME — D) m in E n := B v (x') x [—77, 77] with (x 1 , 0) on its free boundary with outward normal ft , which 
moves with normal velocity 7. Further, S will be an "almost" supersolution near {x 1 ', 0) in the following 
sense: 

7 > \VS\ + VS • V$ - e at (x', 0). 

Proof. Recall that the Barenblatt profiles are given by 

B{x ^ c)= m + rr-Klx-x ^ + 

(t + T) 

where C > and A = ((m — l)n + 2) -1 , K = A/2, and they arc solutions of (PME) m . C and r arc 
parameters that control the free boundary speed and initial support. 

Now we change variables so that ft is colincar with x' , and take xq = . Then it suffices to take 
xq = 0. We start with B(x,t): a Barenblatt solution with initial support Br(0), and initial free 
boundary advancement speed £ where R and £ will be determined later. We fix r(t) = \x — vt, with /j, v 
as yet unspecified. Then we define 

S(x, t) = sup B(y, t) = B((l - r(t)/|a:|)a;, t) in E n , 

yeB r(t) (x) 

where n is for now much smaller than R/2. 

Note that 

S t = B t - r'(t)VB .±-=B t + r'(t)\VB\ =B t - v\VB\ 
\x\ 

Moreover in E Vl since r\ < R/2, then E n is bounded away from the origin and we get l/|x| < 2/R. Thus 
we find that 

dx~ - Wj +^B\0 { 1/R). 

Thus VS 1 = WB + 0(/j,) and since |V5| does not vary fast in E, we can repeat and find AS* = AB + 0(fi). 
Using that B is a (PME) m solution then gives 

S t = B t - u\VB\ = (m- 1)SAS + \WS\ 2 - v\VS\ + 0{p) 

Next let us define 

S(x, t) = S(x + bt, t), where b = V$(x', 0). 

Then S t — S t + V5 • b, and one can conclude that 

S t = (m-l)S'(A5 + A$) + |VS'| 2 + VS'-V$-(m-l)5'A$-^|VS , |+0(M) 

Now in E v , S(x,t) < 2r;sup Br7 |V5| = 0(n), so (m - 1)5A$ = 0(n). Therefore 

S t = (m- 1)S(AS + A*) + |V5| 2 + r'(i)|V5| + V5 • V$ + 0(/*) + 0(77) 

At this point we have to start picking our parameters carefully. First, we can assume that e is small 
enough so that e < inf s n{s>o} l^-SI/S for some small value of n. Then we take 

v = e/3, ^ = 7-6-x' + ^>0 
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Now we take 77, fi small enough so that in E v , 

|O0u) + 0(77)1 <e inf |V5|/3 
Er,n{s>o} 

and we set R = \x'\ — /j = 1 — /i. Now we refine 77 so that 

sup |VS| - inf |V5| < e inf |V5| 
£ r ,n{s>o} £„n{s>o} £„n{s>o} 

Then our choice of v gives us the estimates 

einf Sr ,n{s>o} |V5| 



^|V5| > 



3 

while also 

„|VS| < CBUp ^ n ^>°> |VS| < £ ( inf ^n { ,>Q} IVgj +j = |V5|/3+ e 2 /3 < , inf |V5|/2 

where we used our assumption on e small. Thus we find that 

-e inf IVSI <r'(t)\VS\+O(ri) + O(r(t))<0 

E v n{S>o} 

and so finally 

(m-l)S'(A5 + A$) + |VS , | 2 +VS'-V$-e inf IVSI < 5 t < (m- 1)5(A5 + A$) + |VS| 2 +VS- V$ 

£„n{s>o} 

Then we are done, since it is clear that (x',0) is on the free boundary of 5* and the free boundary has 
initial velocity 7. 

□ 

Appendix B. The support of (PME-D) m solutions have bounded jumps 

Theorem B.l. Suppose u m is a solution to (PME-D) m in R d . Then for K > 0, there exist constants 
i , ma.x,T > only depending on K,d, and $ near x' such that the following holds for any r$ < r max : 
Suppose u m (-,t') = in B ro (x') and u m < K on the parabolic boundary of B2 ro (x') x [t',t' + T]. Then 
we have that u m = in B ro u(x') x [t ', t 1 + T]. 

Proof. 1. We may assume that (x',t') = 0. Now to prove this theorem, first we construct a super- 
solution of (PME-D) m in i?2r o (0) x [0, T] where T is yet to be determined. We start by constructing 
u(r) on i?2r o (0) satisfying u(ro) = K and u(r) = for r < ro/2. We take ro small enough so that 
su Ps 4r (o) |V$(a;) — V$(0)| < 1. Define a = sup B4r ^ A$(x). Then we solve that Aw = —a in the 
annulus B2r (0)\-B ro /2(0). This yields u given by 

_J2 9^1 in a u o 

r d-2 2d -I- >-> a 7= z 

Clnr-ar 2 /4 + Z) N = 2 
We proceed assuming d > 2; the d = 2 case is similar. We choose C and £> so that u(ro) — K, u(ro/2) = 0: 

By taking derivatives it can be seen that u has the largest derivative at ro/2, and we then estimate: 

u'(ro/2) < C(d, JT)/r + C(d, $)dr 2 
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Further, we notice that v! > C^I/t-q -1 ) and so if we take r$ small enough we have u'(r) > 1 when 
r o/2 < r < 4ro. Then this entails that u(r) > u(tq) — K if r € [ro,4ro], and we are done finding r max 
which is the largest value of tq that makes the desired estimates hold. 

2. Now let us define 

u(r,t) :=u{R{t)r) 

where R(t) is a function to be determined with R(0) = 1, 1 < i?(i) < 3/2. Then by construction of u, 

Am = R(t) 2 (Au)(R(t)r) = -aR(tf < -a < 

Further straightforward computation yields that u t > 2|Vu| 2 holds if 

iT(t) u'(R(t)r) 
R(t) 2 r 

To this end, let us choose R(t) = 1/(1 — Li), where 

L:=C(N 7 $,K)/r 2 = - sup u'(r) > 2 sup w l R W r ). 

r re[r /2,4r ] re[r /2,8r /3] r 

3. Lastly we define 

v(x, t) :— u(x + bt, i) 

where b = V$(0,0). We claim that u is a (PME-D) m supersolution in B 2ro (0) x [0, T(d, $, K)] for any 
choice of m. To see this, note that 

u t - 6 • Vv > 2|Vu| 2 > (m - l)u(A« + A$) + |Vu| 2 + |Vu| 2 . 

Now if ro < r max , we have that 

|Vu| 2 - Vu • (V$ - S) > |Vu|(|Vu| - \b - V$| > |Vu|(|Vu| - 1) 

But we know that |Vu| = u'(r) > R(t)v! > u' > 1, so the above quantity is positive. Thus 

v t > (m- l)u(Au + A$) + |Vu| 2 + Vu • V$ 

and since u* > |Vu| 2 + 6- Vu in general this holds at the boundary too. Thus v is a classical free boundary 
supersolution, and so by Lemma 2.6 in |KL| a viscosity supersolution. 



Lastly, T is chosen so that both \b\T < 2r /3 and T < 1/(3L), and lastly so that T < r /(12|6|). The 
first condition ensures that the bounds on u in B Sro / 3 hold for v in B2r a and the second ensures that 
R(t) < 3/2. The last one ensures that if \x\ < r /4, \x + bt\ < \x\ + r /12 < r /3 < r /2R(t) and hence 

v(x, t) = u(R(t)(x + bt), t) = in \x\ < r /4. 

Now by construction u m < v on the parabolic boundary of i?2ro(0) x [0, T(d, <&,K)], so we can apply a 



comparison principle, Theorem 2.25 in KL , to find that 

u m < v in B 2ro (0) x [0,T(d,$,K)} 
Then observing the properties of v, we are done. □ 
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Appendix C. Some prior results on gradient flows 



:S5 



In this part of appendix, we state some results from AGS , concerning the existence and uniqueness 



of the discrete solution p^ n h as defined in ( 4.4 ) , and the convergence as the time step h — > 



The key step leading to these results is the A-convexity of E m [p] for all 1 < m < oo along the 
generalized geodesies. Thus we first digress a little b it to state some definition and results from the 
optimal transport theory (see e.g. section 9.2 in AGS ). Recall that /i £ VilM. ) is regular if /j, £ L p (R d ) 
with some p > 1. 

Definition C.l (generalized geodesies). Let the reference measure p 1 £ 7*2 (U.) be regular. Let p 2 ,p 3 £ 
T'a (K. ) / then we can find two optimal transport maps t and t such that t\p} — p l and W 2 (pi, p l ) = 
L d \t l (x) — x\ dfjr(x) for i = 2,3. The generalized geodesies joining p 2 to p 3 (with base p 1 ) is defined 
as 



/**'- = (^ 3 )#^ where t 2 ^ 3 := (1 - t)t 2 



,2—5-3 



tt 3 , i€[0,l]. 



(C.l) 



Using the notion of generalized geodesies, one can define a notion of semi- convexity (or X-convexity) 
for energy functional on 'Pa(R ) : 

Definition C.2 (A-convexity along generalized geodesies). Given A £ R, a functional E is called X- 
convex along the generalized geodesies if for any pi,p 2 and p 3 satisfying the conditions in Definition 
\C.1\ the following inequality holds 

E[tf^ 3 } < (1 - t)E\jj, 2 ] + tE[fi 3 ] - -t{\ -t) I \t 2 - t 3 | 2 d^i for allO<t<l, 

2 J R d 

where fJ. 2 ^ 3 , t 2 and t 3 are as defined in Definition C.l 



The following Lemma is a direct consequence of AGS Sec 9.3], which says that as long as $ is semi- 
convex, the functional E m would be convex for all 1 < m < oo. Since the case m = oo is not directly 
covered in the book, we provide a short proof below for the sake of completeness. 

Lemma C.3 ( [AGS]). Let $ satisfy (A3) , and let E m : V 2 (R d ) -^ R be as defined as in dO. Then 



E m is X-convex along general geodesies for all 1 < m < oo. 



Proof. Due to (A3), Proposition 9.3.2 of AGS gives the A-convexity of the functional L d p$dx along 
generalized geodesies. 



For a finite m > 1, let S m be given by (4.2). One can directly apply Proposition 9.3.9 in |AGS to 
obtain the convexity of S m along generalized geodesies. Since the sum of two A-convex functionals is 
still A convex, we obtain the A-convexity of E m for any finite m > 1. 



It remains to check that the functional Soo defined in 
To do this, let //, i = 1,2,3 be as given in Definition 



4.3) is also A-convex along generalized geodesies. 



C.l It suffices to show that if ||//||z«> < 1 for 



i = 2, 3, then || A^f - *" 3 || j l°= < 1 for all < t < 1 as well. Note that due to the A-convexity of S m for all 



m > 1, we obtain 

||A^ 3 |k™ < min{||^ 2 || Lm , \\p 3 \\ L ™} for all m > 1, 
and sending to — > oo immediately yields the desired result. 



□ 



Once we have the A-convexity of E m , Lemma 9.2.7 in AGS guarantees that the Assumption 4.0.1 



in AGS is satisfied, which leads to the existence, uniqueness and convergence results in Theorem 4.1 
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